Gene expression, evolution, and the genetics of electrosensing in the smalltooth sawfish, Pristis pectinata.

Sawfishes (Pristidae) are large, highly threatened rays named for their tooth-studded rostrum, which is used for prey sensing and capture. Of all five species, the smalltooth sawfish, Pristis pectinata, has experienced the greatest decline in range, currently found in only ~20% of its historic range. To better understand the genetic underpinnings of these taxonomically and morphologically unique animals, we collected transcriptomic data from several tissue types, mapped them to the recently completed reference genome, and contrasted the patterns observed with comparable data from other elasmobranchs. Evidence of positive selection was detected in 79 genes in P. pectinata, several of which are involved in growth factor/receptor tyrosine kinase signaling and body symmetry and may be related to the unique morphology of sawfishes. Changes in these genes may impact cellular responses to environmental conditions such as temperature, dissolved oxygen, and salinity. Data acquired also allow for examination of the molecular components of P. pectinata electrosensory systems, which are highly developed in sawfishes and have likely been influential in their evolutionary success.

Insights into the nuclear and mitochondrial genome of the Lemon shark Negaprion brevirostris using low-coverage sequencing: Genome size, repetitive elements, mitochondrial genome, and phylogenetic placement.

The Lemon shark Negaprion brevirostris is an important species experiencing conservation issues that is in need of genomic resources. Herein, we conducted a genome survey sequencing in N. brevirostris and determined genome size, explored repetitive elements, assembled and annotated the 45S rRNA DNA operon, and assembled and described in detail the mitochondrial genome. Lastly, the phylogenetic position of N. brevirostris in the family Carcharhinidae was examined using translated protein coding genes. The estimated haploid genome size ranged between 2.29 and 2.58 Gbp using a k-mer analysis, which is slightly below the genome size estimated for other sharks belonging to the family Carcharhinidae. Using a k-mer analysis, approx. 64-71 % of the genome of N. brevirostris was composed of repetitive elements. A relatively large proportion of the 'repeatome' could not be annotated. Taking into account only annotated repetitive elements, Class I - Long Interspersed Nuclear Element (LINE) were the most abundant repetitive elements followed by Class I - Penelope and Satellite DNA. The nuclear ribosomal operon was fully assembled. The AT-rich complete mitochondrial genome was 16,703 bp long and encoded 13 protein coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Negaprion brevirostris is closely related to the genera Carcharhinus, Glyphis and Lamiopsis in the family Carcharinidae. This new genomic resources will aid with the development of conservation plans for this large coastal shark.

Accounting for unobserved population dynamics and aging error in close-kin mark-recapture assessments.

Obtaining robust estimates of population abundance is a central challenge hindering the conservation and management of many threatened and exploited species. Close-kin mark-recapture (CKMR) is a genetics-based approach that has strong potential to improve the monitoring of data-limited species by enabling estimates of abundance, survival, and other parameters for populations that are challenging to assess. However, CKMR models have received limited sensitivity testing under realistic population dynamics and sampling scenarios, impeding the application of the method in population monitoring programs and stock assessments. Here, we use individual-based simulation to examine how unmodeled population dynamics and aging uncertainty affect the accuracy and precision of CKMR parameter estimates under different sampling strategies. We then present adapted models that correct the biases that arise from model misspecification. Our results demonstrate that a simple base-case CKMR model produces robust estimates of population abundance with stable populations that breed annually; however, if a population trend or non-annual breeding dynamics are present, or if year-specific estimates of abundance are desired, a more complex CKMR model must be constructed. In addition, we show that CKMR can generate reliable abundance estimates for adults from a variety of sampling strategies, including juvenile-focused sampling where adults are never directly observed (and aging error is minimal). Finally, we apply a CKMR model that has been adapted for population growth and intermittent breeding to two decades of genetic data from juvenile lemon sharks (Negaprion brevirostris) in Bimini, Bahamas, to demonstrate how application of CKMR to samples drawn solely from juveniles can contribute to monitoring efforts for highly mobile populations. Overall, this study expands our understanding of the biological factors and sampling decisions that cause bias in CKMR models, identifies key areas for future inquiry, and provides recommendations that can aid biologists in planning and implementing an effective CKMR study, particularly for long-lived data-limited species.

Parthenogenesis in an elasmobranch in the presence of conspecific males.

Parthenogenesis has been observed in several elasmobranch species, primarily in public aquaria. The majority of cases of parthenogenesis have occurred either when females were held without males or once a male was removed from a female's habitat. Here we report a second instance of parthenogenesis in a zebra shark female that was housed with conspecific mature males. This study calls into question the conditions under which elasmobranch females undergo parthenogenesis.

Deep ecomorphological and genetic divergence in Steller's Jays (Cyanocitta stelleri, Aves: Corvidae).

The relationship between ecology and morphology is a cornerstone of evolutionary biology, and quantifying variation across environments can shed light on processes that give rise to biodiversity. Three morphotypes of the Steller's Jay (Cyanocitta stelleri) occupy different ecoregions in western North America, which vary in climate and landcover. These morphotypes (Coastal, Interior, Rocky Mountain) differ in size, plumage coloration, and head pattern. We sampled 1080 Steller's Jays from 68 populations (plus 11 outgroups) to address three main questions using data on morphology, plumage, genetics (mtDNA, microsatellites), and ecological niches: (1) How do phenotypic and genetic traits vary within and among populations, morphotypes, and ecoregions? (2) How do population-level differences in Steller's Jays compare with other sister species pairs of North American birds? (3) What can we infer about the population history of Steller's Jays in relation to past climates, paleoecology, and niche evolution? We found substantial morphological, genetic, and ecological differentiation among morphotypes. The greatest genetic divergence separated Coastal and Interior morphotypes from the Rocky Mountain morphotype, which was associated with warmer, drier, and more open habitats. Microsatellites revealed additional structure between Coastal and Interior groups. The deep mtDNA split between Coastal/Interior and Rocky Mountain lineages of Steller's Jay (ND2 ~ 7.8%) is older than most North American avian sister species and dates to approximately 4.3 mya. Interior and Rocky Mountain morphotypes contact across a narrow zone with steep clines in traits and reduced gene flow. The distribution of the three morphotypes coincides with divergent varieties of ponderosa pine and Douglas fir. Species distribution models support multiple glacial refugia for Steller's Jays. Our integrative dataset combined with extensive geographic sampling provides compelling evidence for recognizing at least two species of Steller's Jay.

Genetic population dynamics of the critically endangered scalloped hammerhead shark (Sphyrna lewini) in the Eastern Tropical Pacific.

The scalloped hammerhead shark, Sphyrna lewini, is a Critically Endangered, migratory species known for its tendency to form iconic and visually spectacular large aggregations. Herein, we investigated the population genetic dynamics of the scalloped hammerhead across much of its distribution in the Eastern Tropical Pacific (ETP), ranging from Costa Rica to Ecuador, focusing on young-of-year animals from putative coastal nursery areas and adult females from seasonal aggregations that form in the northern Galápagos Islands. Nuclear microsatellites and partial mitochondrial control region sequences showed little evidence of population structure suggesting that scalloped hammerheads in this ETP region comprise a single genetic stock. Galápagos aggregations of adults were not comprised of related individuals, suggesting that kinship does not play a role in the formation of the repeated, annual gatherings at these remote offshore locations. Despite high levels of fisheries exploitation of this species in the ETP, the adult scalloped hammerheads here showed greater genetic diversity compared with adult conspecifics from other parts of the species' global distribution. A phylogeographic analysis of available, globally sourced, mitochondrial control region sequence data (n = 1818 sequences) revealed that scalloped hammerheads comprise three distinct matrilines corresponding to the three major world ocean basins, highlighting the need for conservation of these evolutionarily unique lineages. This study provides the first view of the genetic properties of a scalloped hammerhead aggregation, and the largest sample size-based investigation of population structure and phylogeography of this species in the ETP to date.

Phylogeography and population genetics of the cryptic bonnethead shark Sphyrna aff. tiburo in Brazil and the Caribbean inferred from mtDNA markers.

Resolving the identity, phylogeny and distribution of cryptic species within species complexes is an essential precursor to management. The bonnethead shark, Sphyrna tiburo, is a small coastal shark distributed in the Western Atlantic from North Carolina (U.S.A.) to southern Brazil. Genetic analyses based on mitochondrial markers revealed that bonnethead sharks comprise a species complex with at least two lineages in the Northwestern Atlantic and the Caribbean (S. tiburo and Sphyrna aff. tiburo, respectively). The phylogeographic and phylogenetic analysis of two mitochondrial markers [control region (mtCR) and cytochrome oxidase I (COI)] showed that bonnethead sharks from southeastern Brazil correspond to S. aff. tiburo, extending the distribution of this cryptic species >5000 km. Bonnethead shark populations are only managed in the U.S.A. and in the 2000s were considered to be regionally extinct or collapsed in southeast Brazil. The results indicate that there is significant genetic differentiation between S. aff. tiburo from Brazil and other populations from the Caribbean (ΦST  = 0.9053, P < 0.000), which means that collapsed populations in the former are unlikely to be replenished from Caribbean immigration. The species identity of bonnethead sharks in the Southwest Atlantic and their relationship to North Atlantic and Caribbean populations still remains unresolved. Taxonomic revision and further sampling are required to reevaluate the status of the bonnethead shark complex through its distribution range.

Artificial insemination and parthenogenesis in the whitespotted bamboo shark Chiloscyllium plagiosum.

Non-lethal methods for semen collection from elasmobranchs to better understand species reproduction has accompanied the development of artificial insemination. Ejaculates (n = 82) collected from whitespotted bamboo sharks Chiloscyllium plagiosum (n = 19) were assessed and cold-stored raw or extended at 4 °C. Females (n = 20) were inseminated with fresh or 24-48 h cold-stored raw or extended semen and paternity of offspring determined with microsatellite markers. Insemination of females with fresh semen (n = 10) resulted in 80 hatchlings and 27.6% fertility. Insemination of females with semen cold-stored 24 h (n = 4) and 48 h (n = 1) semen resulted in 17 hatchlings and fertilization rates of 28.1% and 7.1% respectively. Two females inseminated with fresh or cold-stored semen laid eggs that hatched from fertilization and parthenogenesis within the same clutch. Parthenogenesis rate for inseminated females was 0.71%. Results demonstrate artificial insemination with cold-stored semen can provide a strategy for transport of male genetics nationally and internationally, precluding the need to transport sharks. Production of parthenotes in the same clutch as sexually fertilized eggs highlights the prevalence of parthenogenesis in whitespotted bamboo sharks and poses important considerations for population management.

Annual breeding in a captive smalltooth sawfish, Pristis pectinata.

The critically endangered smalltooth sawfish Pristis pectinata reproduces biennially in central west Florida, U.S.A. Here we demonstrate that smalltooth sawfish are physiologically capable of reproducing annually in a captive environment. The smalltooth sawfish are held in an open system, with abiotic conditions that vary naturally with the surrounding environment in The Bahamas. This suggests wild smalltooth sawfish may also be capable of annual reproduction provided there are adequate prey resources, limited competition and mate availability.

Species composition of the largest shark fin retail-market in mainland China.

Species-specific monitoring through large shark fin market surveys has been a valuable data source to estimate global catches and international shark fin trade dynamics. Hong Kong and Guangzhou, mainland China, are the largest shark fin markets and consumption centers in the world. We used molecular identification protocols on randomly collected processed fin trimmings (n = 2000) and non-parametric species estimators to investigate the species composition of the Guangzhou retail market and compare the species diversity between the Guangzhou and Hong Kong shark fin retail markets. Species diversity was similar between both trade hubs with a small subset of species dominating the composition. The blue shark (Prionace glauca) was the most common species overall followed by the CITES-listed silky shark (Carcharhinus falciformis), scalloped hammerhead shark (Sphyrna lewini), smooth hammerhead shark (S. zygaena) and shortfin mako shark (Isurus oxyrinchus). Our results support previous indications of high connectivity between the shark fin markets of Hong Kong and mainland China and suggest that systematic studies of other fin trade hubs within Mainland China and stronger law-enforcement protocols and capacity building are needed.

Metagenomic 16S rDNA amplicon data on bacterial diversity profiling and its predicted metabolic functions of varillales in Allpahuayo-Mishana National Reserve.

The white-sands forests or varillales of the Peruvian Amazon are characterized by their distinct physical characteristics, patchy distribution, and endemism [1, 2]. Much research has been conducted on the specialized plant and animal communities that inhabit these ecosystems, yet their soil microbiomes have yet to be studied. Here we provide metagenomic 16S rDNA amplicon data of soil microbiomes from three types of varillales in Allpahuayo-Mishana National Reserve near Iquitos, Peru. Composite soil samples were collected from very low varillal, high-dry varillal, and high-wet varillal. Purified metagenomic DNA was used to prepare and sequence 16S rDNA metagenomic libraries on the Illumina MiqSeq platform. Raw paired-endsequences were analyzed using the Metagenomics RAST server (MG-RAST) and Parallel-Meta3 software and revealed the existence of a high percentage of undiscovered sequences, potentially indicating specialized bacterial communities in these forests. Also, were predicted several metabolic functions in this dataset. The raw sequence data in fastq format is available in the public repository Discover Mendeley Data (https://data.mendeley.com/datasets/syktzxcnp6/2). Also, is available at NCBI's Sequence Read Archive (SRA) with accession numbers SRX7891206 (very low varillal), SRX7891207 (high-dry varillal), and SRX7891208 (high-wet varillal).

Microsatellite primer development in elasmobranchs using next generation sequencing of enriched libraries.

Microsatellites are useful in studies of population genetics, sibship, and parentage. Here, we screened for microsatellites from multiple elasmobranch genomic libraries using an enrichment protocol followed by sequencing on an Illumina platform. We concurrently screened five and then nine genomes and describe the number of potential loci from each respective round of sequencing. To validate the efficacy of the protocol, we developed and tested primers for the pelagic thresher shark, Alopias pelagicus. The method described here is a cost-effective protocol to increase the pool of potential useful loci and allows the concurrent screening of multiple libraries.

No apparent negative tagging effects after 13 years at liberty for lemon shark, Negaprion brevirostris implanted with acoustic transmitter.

An intact and uncompromised internal acoustic transmitter was non-lethally recovered from a lemon shark Negaprion brevirostris, after 13 years at liberty. The shark, first tagged at an estimated age of 2 years old near South Bimini, Bahamas in 2004, was recaptured in 2017 with a total length of 264 cm. The tagged shark displayed typical growth rate, pregnancy, natal homing and pupping behaviour of other individuals in this population. This observation provides important evidence regarding the effects from long-term retention of implanted acoustic transmitters in a carcharhinid shark.

Species composition of the international shark fin trade assessed through a retail-market survey in Hong Kong.

The shark fin trade is a major driver of shark exploitation in fisheries all over the world, most of which are not managed on a species-specific basis. Species-specific trade information highlights taxa of particular concern and can be used to assess the efficacy of management measures and anticipate emerging threats. The species composition of the Hong Kong Special Administrative Region of China, one of the world's largest fin trading hubs, was partially assessed in 1999-2001. We randomly selected and genetically identified fin trimmings (n = 4800), produced during fin processing, from the retail market of Hong Kong in 2014-2015 to assess contemporary species composition of the fin trade. We used nonparametric species estimators to determine that at least 76 species of sharks, batoids, and chimaeras supplied the fin trade and a Bayesian model to determine their relative proportion in the market. The diversity of traded species suggests species substitution could mask depletion of vulnerable species; one-third of identified species are threatened with extinction. The Bayesian model suggested that 8 species each comprised >1% of the fin trimmings (34.1-64.2% for blue [Prionace glauca], 0.2-1.2% for bull [Carcharhinus leucas] and shortfin mako [Isurus oxyrinchus]); thus, trade was skewed to a few globally distributed species. Several other coastal sharks, batoids, and chimaeras are in the trade but poorly managed. Fewer than 10 of the species we modeled have sustainably managed fisheries anywhere in their range, and the most common species in trade, the blue shark, was not among them. Our study and approach serve as a baseline to track changes in composition of species in the fin trade over time to better understand patterns of exploitation and assess the effects of emerging management actions for these animals.

Risky business for a juvenile marine predator? Testing the influence of foraging strategies on size and growth rate under natural conditions.

Mechanisms driving selection of body size and growth rate in wild marine vertebrates are poorly understood, thus limiting knowledge of their fitness costs at ecological, physiological and genetic scales. Here, we indirectly tested whether selection for size-related traits of juvenile sharks that inhabit a nursery hosting two dichotomous habitats, protected mangroves (low predation risk) and exposed seagrass beds (high predation risk), is influenced by their foraging behaviour. Juvenile sharks displayed a continuum of foraging strategies between mangrove and seagrass areas, with some individuals preferentially feeding in one habitat over another. Foraging habitat was correlated with growth rate, whereby slower growing, smaller individuals fed predominantly in sheltered mangroves, whereas larger, faster growing animals fed over exposed seagrass. Concomitantly, tracked juveniles undertook variable movement behaviours across both the low and high predation risk habitat. These data provide supporting evidence for the hypothesis that directional selection favouring smaller size and slower growth rate, both heritable traits in this shark population, may be driven by variability in foraging behaviour and predation risk. Such evolutionary pathways may be critical to adaptation within predator-driven marine ecosystems.

Global population genetic dynamics of a highly migratory, apex predator shark.

Knowledge of genetic connectivity dynamics in the world's large-bodied, highly migratory, apex predator sharks across their global ranges is limited. One such species, the tiger shark (Galeocerdo cuvier), occurs worldwide in warm temperate and tropical waters, uses remarkably diverse habitats (nearshore to pelagic) and possesses a generalist diet that can structure marine ecosystems through top-down processes. We investigated the phylogeography and the global population structure of this exploited, phylogenetically enigmatic shark by using 10 nuclear microsatellites (n = 380) and sequences from the mitochondrial control region (CR, n = 340) and cytochrome oxidase I gene (n = 100). All three marker classes showed the genetic differentiation between tiger sharks from the western Atlantic and Indo-Pacific ocean basins (microsatellite FST  > 0.129; CR ΦST  > 0.497), the presence of North vs. southwestern Atlantic differentiation and the isolation of tiger sharks sampled from Hawaii from other surveyed locations. Furthermore, mitochondrial DNA revealed high levels of intraocean basin matrilineal population structure, suggesting female philopatry and sex-biased gene flow. Coalescent- and genetic distance-based estimates of divergence from CR sequences were largely congruent (dcorr  = 0.0015-0.0050), indicating a separation of Indo-Pacific and western Atlantic tiger sharks <1 million years ago. Mitochondrial haplotype relationships suggested that the western South Atlantic Ocean was likely a historical connection for interocean basin linkages via the dispersal around South Africa. Together, the results reveal unexpectedly high levels of population structure in a highly migratory, behaviourally generalist, cosmopolitan ocean predator, calling for management and conservation on smaller-than-anticipated spatial scales.

Surgeons and suture zones: Hybridization among four surgeonfish species in the Indo-Pacific with variable evolutionary outcomes.

Closely related species can provide valuable insights into evolutionary processes through comparison of their ecology, geographic distribution and the history recorded in their genomes. In the Indo-Pacific, many reef fishes are divided into sister species that come into secondary contact at biogeographic borders, most prominently where Indian Ocean and Pacific Ocean faunas meet. It is unclear whether hybridization in this contact zone represents incomplete speciation, secondary contact, an evolutionary dead-end (for hybrids) or some combination of the above. To address these issues, we conducted comprehensive surveys of two widely-distributed surgeonfish species, Acanthurus leucosternon (N=141) and A. nigricans (N=412), with mtDNA cytochrome b sequences and ten microsatellite loci. These surgeonfishes are found primarily in the Indian and Pacific Oceans, respectively, but overlap at the Christmas and Cocos-Keeling Islands hybrid zone in the eastern Indian Ocean. We also sampled the two other Pacific members of this species complex, A. achilles (N=54) and A. japonicus (N=49), which are known to hybridize with A. nigricans where their ranges overlap. Our results indicate separation between the four species that range from the recent Pleistocene to late Pliocene (235,000-2.25million years ago). The Pacific A. achilles is the most divergent (and possibly ancestral) species with mtDNA dcorr≈0.04, whereas the other two Pacific species (A. japonicus and A. nigricans) are distinguishable only at a population or subspecies level (ΦST=0.6533, P<0.001). Little population structure was observed within species, with evidence of recent population expansion across all four geographic ranges. We detected sharing of mtDNA haplotypes between species and extensive hybridization based on microsatellites, consistent with later generation hybrids but also the effects of allele homoplasy. Despite extensive introgression, 98% of specimens had concordance between mtDNA lineage and species identification based on external morphology, indicating that species integrity may not be eroding. The A. nigricans complex demonstrates a range of outcomes from incomplete speciation to secondary contact to decreasing hybridization with increasing evolutionary depth.

Facultative parthenogenesis in a critically endangered wild vertebrate.

Facultative parthenogenesis - the ability of sexually reproducing species to sometimes produce offspring asexually - is known from a wide range of ordinarily sexually reproducing vertebrates in captivity, including some birds, reptiles and sharks [1-3]. Despite this, free-living parthenogens have never been observed in any of these taxa in the wild, although two free-living snakes were recently discovered each gestating a single parthenogen - one copperhead (Agkistrodon contortrix) and one cottonmouth (Agkistrodon piscivorus) [1]. Vertebrate parthenogens are characterized as being of the homogametic sex (e.g., females in sharks, males in birds) and by having elevated homozygosity compared to their mother [1-3], which may reduce their viability [4]. Although it is unknown if either of the parthenogenetic snakes would have been carried to term or survived in the wild, facultative parthenogenesis might have adaptive significance [1]. If this is true, it is reasonable to hypothesize that parthenogenesis would be found most often at low population density, when females risk reproductive failure because finding mates is difficult [5]. Here, we document the first examples of viable parthenogens living in a normally sexually reproducing wild vertebrate, the smalltooth sawfish (Pristis pectinata). We also provide a simple approach to screen any microsatellite DNA database for parthenogens, which will enable hypothesis-driven research on the significance of vertebrate parthenogenesis in the wild.

Isolation and development of microsatellite loci in an African Woodpecker (Campethera nivosa) using polymerase chain reaction and DNA sequencing.

BACKGROUND: The Buff-spotted Woodpecker (Campethera nivosa) is a resident bird species that is distributed in lowland rainforest habitats from western to eastern Africa. We developed species-specific microsatellite markers to examine the population genetics of this species. FINDINGS: Twenty-one microsatellite loci were isolated from C. nivosa. Of these, 15 were found to amplify consistently. These loci were then tested for variability in 15 individuals from different lowland forest localities. The number of alleles ranged from 3 to 13 per locus, with observed and expected heterozygosity ranging from 0.100 to 0.917 and 0.485 to 0.901, respectively. Four loci exhibited significant heterozygote deficiency while one had an excess of heterozygotes. None of the loci exhibited linkage disequilibrium. CONCLUSION: These polymorphic microsatellite markers will be used to study genetic variability in populations of C. nivosa across either sides of the Congo River to evaluate the effect of the river as a barrier to gene flow.

A quantitative genetic approach to assess the evolutionary potential of a coastal marine fish to ocean acidification.

Assessing the potential of marine organisms to adapt genetically to increasing oceanic CO2 levels requires proxies such as heritability of fitness-related traits under ocean acidification (OA). We applied a quantitative genetic method to derive the first heritability estimate of survival under elevated CO2 conditions in a metazoan. Specifically, we reared offspring, selected from a wild coastal fish population (Atlantic silverside, Menidia menidia), at high CO2 conditions (∼2300 µatm) from fertilization to 15 days posthatch, which significantly reduced survival compared to controls. Perished and surviving offspring were quantitatively sampled and genotyped along with their parents, using eight polymorphic microsatellite loci, to reconstruct a parent-offspring pedigree and estimate variance components. Genetically related individuals were phenotypically more similar (i.e., survived similarly long at elevated CO2 conditions) than unrelated individuals, which translated into a significantly nonzero heritability (0.20 ± 0.07). The contribution of maternal effects was surprisingly small (0.05 ± 0.04) and nonsignificant. Survival among replicates was positively correlated with genetic diversity, particularly with observed heterozygosity. We conclude that early life survival of M. menidia under high CO2 levels has a significant additive genetic component that could elicit an evolutionary response to OA, depending on the strength and direction of future selection.